Aharonov-Bohm cages in two-dimensional structures 
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We present an extreme localization mechanism induced by a magnetic field for tight-binding elec- 
trons in two-dimensional structures. This spectacular phenomenon is investigated for a large class 
of tilings (periodic, quasiperiodic, or random). We are led to introduce the Aharonov-Bohm cages 
defined as the set of sites eventually visited by a wavepacket that can, for particular values of the 
magnetic flux, be bounded. We finally discuss the quantum dynamics which exhibits an original 
pulsating behaviour. 

PACS numbers: 7100, 7210, 7335. 
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Over the past two decades, the behaviour of electrons 
in two-dimensional structures under a magnetic field has 
been of special interest in condensed matter physics. In 
situations where electron-electron interactions are domi- 
nant (low carrier density heterostructures for instance), a 
strong magnetic field has been shown to induce new col- 
lective states such as fractional quantum Hall states 
or Wigner crystals . But even when electrons may be 
considered as weakly interacting, fascinating phenomena 
may take place in the presence of a spatially modulated 
potential besides the external magnetic field. In the case 
of a random potential, one obtains the integer quantum 
Hall effect || which is characterized by a complex one- 
particle spectrum of localized states and a discrete set of 
eigenenergies associated to extended states If the 
modulated potential corresponds to a periodic lattice, 
a subtle competition between this periodicity and the 
area scale imposed by the magnetic field induces remark- 
able features for the one-particle spectrum || , illustrated 
by the well-known Hofstadter's butterfly-like pattern ||. 
The physics of this magnetically induced frustration has 
also generated many experiments on submicronic struc- 
tures such as superconducting networks 0. They have 
in turn triggered investigations of several different two- 
dimensional structures, like the triangular lattic e [p] , the 
honeycomb tiling Jj|, the quasiperiodic tilings |l(fl , and 
even fractal networks |tlj (l2) . 

In this letter, we discuss a new effect of extreme local- 
ization in a large class of rhombus tilings. For particu- 
lar values of the magnetic field, the set of sites visited 
by an initially localized wavepacket can be bounded due 
to Aharonov-Bohm (A.B.) destructive interferences [ p^[ . 
We are thus led to introduce the notion of A.B. cages. 
Our aim is to analyze the relationship between this phe- 
nomenon, the tiling geometry, the quantum dynamics, 
and the spectrum. We therefore pay attention to a peri- 
odic tiling with hexagonal symmetry pi], here denoted 
T 3 (see figure ^). We show that its spectrum is related to 
that of the triangular lattice through a simple analytical 



transformation, which allows to compute the eigenspec- 
trum for rational values of the magnetic flux. When the 
magnetic flux per rhombic plaquette equals half the flux 
quantum, one faces a rather exotic effect for an infinite 
periodic tiling, since the density of states degenerates into 
three 5-peaks and the electron motion becomes bounded. 
We also present another example of an infinite periodic 
tiling denoted T^, which, while having a more complex 
geometry, displays the same properties as 7%. 
Other structures, like quasiperiodic or random tilings, 
can contain some bounded cages due to these particu- 
lar star configurations, but generically, one does not ex- 
pect a highly degenerated discrete spectrum as for T 3 and 
T4. In both cases, cages induce interesting instances of 
wavepacket time evolution with well characterized puls- 
ing states. 

Let us first examine T3, a bipartite periodic hexago- 
nal structure with 3 sites per unit cell (figure [l]), one 
6-fold coordinated (denoted A) and two 3-fold coordi- 
nated (denoted B and C). We consider an tight-binding 
hamiltonian defined by: 



(i) 



<i,j> 



where tij = 1 if i and j are nearest neighbours, oth- 
erwise, and \i > is a localized orbital on site i. In the 
presence of a magnetic field H, the hopping terms are 
modified by phase factors involving the vector potential 
A |15|. Let jij be the phase factor between sites i and j: 
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where 0o = hc/e is the flux quantum. We first con- 
sider a uniform magnetic field H = Hz in the Landau 
gauge: A = H(0, x, 0). As shown below, the whole spec- 
trum only depends on the reduced flux / = <p/4>o where 
4> = Ha 2 ^/3/2 is the magnetic flux through an elemen- 
tary rhombus. Taking into account the translational in- 
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variance in the y direction due to the gauge choice, the 
eigenfunctions write: 



<p(x,y) = ip{x)e lk y y , k y £ [o,27r/aV3 



(3) 



which allows to map the spectrum onto the solution of 
an effective one-dimensional problem: 



(e z -6)yj m = 2 cos ^ 



2r(is I y ( m + 5 ) + " ) ' '»'■-> 
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+2 cos (37m + 2k) ip„ 



(4) 



where tp m = ^Pa(x) for x = ima/2 (m £ Z), is the wave 
function at the A-type sites, K — k y \/3/2, 7 = 2nf, and 
e is the energy for the two-dimensional original problem. 
Note that this mapping is only valid for e ^ (the e = 
case is briefly adressed below). For the triangular lattice, 
a similar equation has been derived by Claro and Wannier 



£T^m = 2 cos ( 7 T ( m + - ) + K ) 1p m + I 



+2 cos |^7 T I m - -J + /e ) Vm-i 
+2 cos(27xm + 2n)ip m 
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with similar notations. So, the 73 spectrum can be 
obtained from the triangular lattice one by choosing 
It = 37/2; the eigenenergies of the two tilings being 
related by: 



6 = 2 cos 



(i) 



(6) 



This extends a classical result on bipartite structures 
to the case of a uniform magnetic field. Since the 
tiling is made up of rhombi, the spectrum is symmet- 
ric (e <-> — e). Moreover, eq. (||) displays a translation 
symmetry / — > / + n (n € N) and a reflection invariance 
about half-integer values of /. Thanks to these symme- 
tries, we can limit our analysis to < / < 1/2. For ratio- 
nal values of / = p/q (p,q mutually prime), the system 
(||) becomes closed after a translation by q periods and 
the spectrum is made up of 2q bands. Note that, when 
q = 3q' (q' £ N) , this period is reduced by a factor 3 (and 
the spectrum has only 2q' bands). An interesting exam- 
ple occurs for / = 1/3 where the spectrum displays two 
symmetric bands extending from ±3 to 0, and is there- 
fore gapless. In addition, for any value of the magnetic 
field, a close inspection of the secular equations system 
for each type of site shows that the spectrum contains a 
highly degenerated eigenenergy (flat band) at e = with 
weight 1/3. Figure ^| shows the 7jj spectrum support ver- 
sus /. The most remarkable point in this context, is that 



1/2, the spectrum collapses into three eigenval- 
and e ± = ±V6 (see eq. 
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One way to analyze the link between the spectrum and 
the quantum dynamics is to study the spectral properties 
at a local level. The magnetic field being uniform, the 
whole spectrum is recovered from the three local density 
of states (LDOS); furthermore, by symmetry, the two 3- 
fold coordinated sites (B and C) are equivalent, so that 
we just need to study two different types of local environ- 
ment. It is then suitable to shift from the above Landau 
gauge to an equivalent cylindrical symmetric gauge given 
byA=(H/2)(-y,x,0). 

We proceed to an analytic Lanczos tridiagonalization of 
local clusters, along the recursion method algorithm [^6| . 
An effective semi-infinite chain is generated, with normal- 
ized orbitals \n) and hopping terms (a n ,b n ), according to 
the three-term recurrence relation: 



b n+ i \n + 1) = (H - a n ) \n) - b n \n - 1) 



(7) 



which allows to evaluate the LDOS at a given site (or any 
linear combination of site orbitals) . Since TL is purely off- 
diagonal and T3 is bipartite (and therefore only contains 
even loops), the diagonal a n terms vanish (i.e. the odd 
moments). Notice that in this case, the orbital |n) en- 
compasses all the sites reached from |0) at least after n 
jumps, and can be interpreted in terms of shell orbitals. 
If one of the b n vanishes, the effective chain becomes fi- 
nite and the LDOS is discrete. It also implies that a 
wavepacket initially localized in the tiling at the origin 
site will only spread over a finite set of sites belonging to 
the mutually coupled shells. 

We first choose the initial orbital |0) at a 3-fold coordi- 
nated site and apply the recursion procedure. One then 
easily obtains: b\ = \/3, b\ = D 2 /VS and 63 = D 3 /D 2 
where: 



D 2 = (12cos 2 (7r/)+9) 1/2 

D 3 = (24cos 2 (tt/) + 12(cos(tt/) + cos(3tt/)) 2 ) 1/2 . (8) 

It is clear that b^ = for / = 1/2, and that the LDOS 
reduces to the following three eigenvalues 0, ±y6- A 
similar derivation for the 6-fold coordined site leads to: 
b\ = VE,b% = 2cos(tt/). For / = 1/2, b\ also vanishes 
and the LDOS reduces to ±\/6 (see figure ||). 
We now proceed to describe the quantum dynamics on 
73. Owing to the superposition principle, we only con- 
sider initial distributions strictly localized on each type 
site, given that any wavepacket behaviour can be de- 
scribed in terms of these individual evolutions. For a 
generic value of /, the initially localized electron wave- 
function spreads with time and asymptotically vanishes 
at the origin, the time dependance of this phenomenon 
being determined by the nature of the eigenspectrum. 
For / = 1/2, the quantum dynamics can be derived in 
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closed forms and it displays a very different behaviour 
since the system remains localized inside a restricted re- 
gion. Indeed, an initial distribution localized on a A 
site is completely "trapped" inside its surrounding star 
(Kj = 0) so that, at any time, there will only be nonvan- 
ishing components of the wavefunction at the origin and 
on the first shell. This oscillating motion also occurs for a 
wavepacket initially localized on a B (or C) site (63 = 0), 
but is here confined inside a larger neighborhood (see fig- 
ure ^|). Consequently, for / = 1/2, any wavepacket on 73 
has a limited extent. 

This extreme localization in a magnetic field can be sim- 
ply understood in terms of Aharonov-Bohm interferences 
around the rhombi; we therefore call the above confined 
regions A.B. cages. More precisely, for a given tiling em- 
bedded in a magnetic field H, we associate to any site 
an A.B. cage Cf defined as the set of sites reached 
by a wavepacket initially localized on Note that 
is generically infinite, but specific values of the magnetic 
field (/ = 1/2 for Tj) can leave it bounded. In this case, 
the cage can also be characterized by its boundary sites, 
where completely destructive interferences disconnect it 
from the rest of the tiling. For a given structure, it is 
sufficient to characterize the inequivalent site cages, and 
several cases can occur: 

• All cages are unbounded for any /. 

Ex: Square lattice, triangular lattice, honeycomb. 

• Some cages are bounded for particular values of /. 
Ex: Star patterns in Penrose tilings and octagonal 
tilings (see below). 

• For the same values of /, all cages are simultane- 
ously bounded. 

The simplest example of the latter case is encountered 
on 7-j for / = 1/2, but it can also occur for other tilings. 
Let us pay attention to 74 represented figure in ^ which 
is an approximant with 14 sites per unit cell of the octag- 
onal quasiperiodic tiling. This tiling contains three types 
of faces (large square, small square, and parallelogram) , 
but only two different areas. 

The previous recursion analysis can be done on the five 
different types of site (A, B, C, D, E) and one can prove 
that for / = 1/2, all the site cages are bounded (in this 
case / refers to the magnetic flux per small area). Thus, 
the eigenspectrum is also discrete and the quantum dy- 
namics is similar to that discussed on 73 (bounded exten- 
sion of any wavepacket, pulsating states, ...). It is also 
possible to continuously deform 74, keeping the same con- 
nectivity, until all the edges have the same length, and 
the star around A becomes 8-fold symmetric. But any 
metrical transformation changes the nature of the spec- 
trum and the dynamics, so that only some site cages 
remains bounded (C? for / = 1/2). In addition, when 
all the edge lengths are equal, the areas are incommensu- 
rate and the energy spectrum is no longer periodic with /. 



This loss of periodicity does always occur when two tile 
areas are incommensurate as it is the case in quasiperi- 
odic tilings. Indeed, consider the closed paths in the 
tiling contributing, for instance, to the fourth moment of 
the LDOS at a generic site. Among them, only the non 
self- retracing ones are field dependent, and lead to cosine 
terms whose arguments depend on the circumnavigated 
rhombi areas. In those tilings, the spectrum becomes 
then quasiperiodic with / as remarked in |Q [ pT[ . 
Concerning quasiperiodic tilings, some local configura- 
tions display interesting characteristics in the presence 
of a magnetic field. As for 73 and 74, star-like patterns 
can be encountered, providing canonical example of A.B. 
cages. An easy way to understand the properties of such 
an environment is to use the symmetry group of the clus- 
ter. For convenience, we only consider p-fold symmetric 
stars (2p + 1 sites) (see figure || (p=6)), whose first shell 
sites are 3-fold coordinated and where all the edges have 
the same length. In the absence of a magnetic flux, those 
stars remains unchanged under C pv transformations, but 
the magnetic field breaks the reflection invariance. Fi- 
nally, one can only deal with the C p rotation group. We 
denote the p one-dimensional irreductible representations 
by Ti (i = 0,p — 1). If we consider the symmetrical rep- 
resentation To, one can show that the LDOS at the star 
center reduces to e ± = ± v /p, as soon as / = 1/2. This re- 
sult is nothing but a generalization of the above recursion 
study on the Tj A site, to any p-order symmetry. Con- 
sequently, the dynamical properties are also comparable 
and one faces a simple example of A.B. cage that appears 
in quasiperiodic rhombus tilings (Penrose (p = 5, 10), oc- 
tagonal (p = 8)). A close investigation of the other p— 1 
representations shows that when f = m/p (m = 1, p — 1), 
one can identify confined eigenstates associated to e = 0, 
similar to those existing for / = 0. Along the same line, 
one can prove that there exist strictly localized states in 
the clusters studied by Kohmoto et al. j^J, for any val- 
ues of the magnetic field. Some of these states persist 
in the infinite tiling and have actually been numerically 
observed by Schwabe et al. [[[o| . 

It is also important to clarify the effect of disorder on 
this field induced localization. Qucerly, one expects the 
randomness to alter and even destroy the phase matching 
essential for this localization effect. A standard way to 
disorder a system is to randomly modulate the hopping 
terms without changing the tiling configuration; but a 
different approach consist in flipping the rhombic tiles as 
in the Random Tiling Model [19]. In T3, a finite num- 
ber of such flips should neither affect the discrete nature 
of the spectrum nor the confined diffusion process. How- 
ever, a finite density of flips may change those properties. 

Finally, a natural question arises about the possibility 
to experimentally observe this very strong localization 
phenomenon. We should first remark that A.B. cages 
are sensitive to various kinds of disorder which are likely 
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to occur in real systems: fluctuations in the tile areas and 
in the transmission matrix along the edges are strongly 
relevant, but unlikely to fully destroy the localized na- 
ture of the eigenstates. We may just conjecture a strange 
non-monotonous variation of the localization length as a 
function of disorder strength. 

A certainly interesting experiment would be to investi- 
gate the magnetic vortex lattice in a superconducting 
wire network with the geometry. Indeed, the unusual 
nature of the linearized Ginzburg-Landau equation eigen- 
states is to be reflected on the behavior of the solutions 
in the presence of the non-linear term. A crucial issue 
would be to determine whether there is a unique optimal 
lattice structure, or a glassy type of vortex-solid, or even 
a disordered vortex-liquid. These various possibilities are 
still debated in simpler geometries f2(i|| . We believe that 
the systems presented here are good candidates for mul- 
tiple ground state configurations, but the nature of the 
energy barriers connecting them is certainly a challenging 
open problem. 
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Figure Captions 

1. A piece of T3 lattice embedded in a perpendicular 
magnetic field H; a is the rhombus edge length. 

2. T3 butterfly-like energy spectrum versus /. 

3. Effective semi-infinite chain (recursion framework), 
associated to % local environments (6-fold site 
(left), 3-fold site (right)). 

4. A piece of the T4 periodic tiling. The five different 
types of sites are depicted 

(A, B, C, D, E). 
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